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Abstract. We undertake a critical evaluation of recent 
observational information on VL m and £!a in order to iden- 
tify possible sources of systematic errors and effects of 
simplified statistical analyses. We combine observations 
for which the results have been published in the form 
of likelihood contours in the f2 m ,fiA plane. We approx- 
imate the contours by fifth order polynomials, and we 
then use the maximum likelihood method to obtain joint 
likelihood contours for the combined data. In the choice 
of statistical merits we aim at minimum loss of infor- 
mation rather than at minimum variance. We find that 
n = ti m + Qa = 0.99 ± 0.04 ± 0.03, where the first error 
is mainly statistical and the second error is systematical. 
In a flat Universe we find n^ at = 0.31 ± 0.04 ± 0.04. 

Key words: Cosmology: observations - Methods: statis- 
tical 



1. INTRODUCTION 

Our knowledge of the dynamical parameters of the 
Universe describing the cosmic expansion has improved 
rapidly over the last few years, starting with the epochal 
discovery of the large scale anisotropies of the CMB by 
COBE-DMR (Smoot et al. 1992), followed by the dra- 
matic supernova la observations by the High-z Supernova 
Search Team (Riess et al. 1998) and the Supernova Cos- 
mology Project (Perlmutter et al. 1998, 1999), and most 
recently by the measurements of the first acoustic peak 
in the CMB power spectrum in the first results from the 
balloon flights BOOMERANG (de Bernardis et al. 2000) 
and MAXIMA (Balbi et al. 2000, Hanany et al. 2000). 

The list of other recent observations is very long, even 
if one restricts oneself to those having information on both 
the mass density parameter f2 m and the density parameter 



17a of vacuum energy. Recall that £!a is related to the 
cosmological constant A by 



(1) 



Send offprint requests to: S. 
hsrashid@pcu.helsinki.fi 



M. Harun-or-Rashid, e-mail: 



Lineweaver (1998) and Tegmark (1999) have summarized 
and analyzed some 20 more observations of the CMB 
anisotropies (cf. their reference lists). Determinations of 
f2 m and £l\ have been reported from observations on the 
gas fraction in X-ray clusters (Evrard 1997), on X-ray clus- 
ter evolution (Bahcall & Fan 1998, Eke et al. 1998), on the 
cluster mass function and the Lya forest (Weinberg et 
al. 1998), on gravitational lensing (Chiba & Yoshi 1998, 
Helbig 2000, Im et al. 1997), on the Sunyaev-Zel'dovich 
effect (Birkinshaw 1999, Carlstrom et al. 1999), on clas- 
sical double radio sources (Guerra et al. 2000), on galaxy 
peculiar velocities (Zehavi & Dekel 1999), on the evolu- 
tion of galaxies and star creation versus the evolution of 
galaxy luminosity densities (Totani 1997). The large scale 
structure and its power spectrum has been studied in the 
SSRS2 and CfA2 galaxy surveys (da Costa et al., 1994), 
in the Las Campanas Redshift Survey (Schectman et al., 
1996), in the Abell-ACO cluster survey (Retzlaff et al., 
1998), and in the 2dF Galaxy Redshift Survey (Peacock 
et al., 2001). 

Some of the above information has already been used 
to constrain f2 m and 17a, some of it could in principle 
be used that way, but has not been presented in a form 
readily useful to an analyst outside the observer teams. 
It must also be said that much is statistically weak, the 
analyses being simplified and the discussions of possible 
systematic errors absent. 

Nevertheless, the list of large combined data analy- 
ses since 1999 is already long. Lineweaver (1999) com- 
bined the SN la data with CMB data, X-ray cluster 
data, cluster evolution data and double radio sources. 
Le Dour et al. (2000) analyzed only CMB data, whereas 
Tegmark et al. (2000a) combined CMB data with IRAS 
LSS data. Tegmark & Zaldarriaga (2000b, 2000c) and Hu 
et al. (2000) combined BOOMERANG and MAXIMA 
data, Melchiorri et al. (2000) combined BOOMERANG 
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and COBE data. Bridle et al. (2000) combined the CMB 
data with galaxy peculiar velocities and with the SN la 
data. The BOOMERANG, MAXIMA and COBE data 
have been combined with LSS and SN la data by Jaffe 
et al. (2000) and Bond et al. (2000), and combined with 
a different set of LSS data by Novosyadlyj et al. (2000a) 
and Durrer & Novosyadlyj (2000). In a sequence of papers 
we (Roos & Harun-or-Rashid 1998, 1999, 2000) combined 
much of the data quoted above having published an error 
on fl m and using simple x 2 analysis. This of course 
implied believing in the errors and treating them as Gaus- 
sian. 

The conclusion of all these partially overlapping anal- 
yses is that the Universe is consistent with being flat, 
f2o = TO + is near unity, and fl\/fl m is near 2. But 
the analyses differ in methods, in the treatment of errors 
and in their attention to possible systematic errors, so the 
results differ in the precision of these conclusions. 

We undertake here yet another combined study us- 
ing the maximum likelihood method, and paying special 
attention to statistical arguments. In Sec. 2 we discuss 
statistical methods in general and present our method of 
analysis. In Sec. 3 we discuss the data chosen for our anal- 
ysis, all of which have been published graphically as like- 
lihood contours in the fl m ,ClA plane. In Sec. 4 we give 
our estimates for the parameters, and in Sec. 5 we discuss 
some related analyses. 

2. STATISTICAL METHODS 

The observational data generally contain information on a 
large set of parameters. The information about each sin- 
gle parameter is then obtained by fixing some parame- 
ters at known input values, and marginalizing over oth- 
ers. Note that it is quite misleading to report the values 
of each parameter in turn, always carrying out uncondi- 
tional marginalizing over all the others, because then the 
same information has been reused many times. Already 
when marginalizing to obtain the value of the first param- 
eter one has used all the information available. 

There is only one remedy to this: if one is mainly in- 
terested in the values of a small subset of parameters, 
in our case two, one should describe their joint pdf by 
confidence contours (or surfaces or hypersurfaces) in the 
space of those parameters, marginalizing over less inter- 
esting ancillary parameters. The confidence range of the 
second parameter is then conditional on the range of the 
first parameter, and so on. In practice all the conditional 
confidence ranges are determined by the size of the or- 
thogonal box circumscribed around the two- (or higher-) 
dimensional confidence contour. 

In the present analysis we are only interested in the 
values of Q m and Qa, therefore we only use data for which 
the marginalization over ancillary parameters has already 
been carried out. Note that thereby we do use the full 
information of each observation. 



In several observations it has been noted that some pa- 
rameters are strongly correlated. If one marginalizes over 
one of a pair of correlated parameters, the likelihood func- 
tion of the other one becomes quite broad. This is an effect 
we do not try to avoid, because it implies including one 
type of systematic error. 

Let us make a few comments of statistical nature com- 
paring the least squares or \ 2 method with the maximum 
likelihood method. The advantage of the least squares 
method is its simplicity, it is unbiased, and a goodness-of- 
fit value can be obtained by comparing the least squares 
sum with the number of degrees of freedom. The disad- 
vantage is that it requires the pdf of the input data to 
be symmetric, preferably normal, but even then there is 
no guarantee that the final estimate will be normally dis- 
tributed, except asymptotically. To form the least squares 
of very conflicting data is statistically meaningless. 

In contrast, the log- likelihood functions of any data 
can be added, and if there are conflicts, due for instance 
to systematic errors, they will show up as several dips. No 
goodness-of-fit value can be obtained thereby one does 
not risk statistically meaningless statements) , but relative 
confidence levels can be defined. There are no restrictions 
to the symmetry or normality of the pdf of the input data. 
Asymptotically the maximum likelihood estimator attains 
normality faster than the least squares estimator. 

When one compromises between different statistical 
merits, one may conclude that it is more desirable to 
achieve minimum loss of information than minimum vari- 
ance. In the first case one wants to make sure that no other 
single number could contain more information about the 
parameter of interest than the estimate chosen. In the sec- 
ond case one feels that the smaller the variance, the more 
certain one is that the estimate is near the true value. If 
one opts for minimum loss of information, the maximum 
likelihood estimator is optimal in the asymptotic limit. 

A problem that occurs in some of the data we use is 
that the pdf extends into an unphysical region, specifi- 
cally the region Cl m < 0. Even the region Q\ < might 
be considered unphysical in this context. To simply ignore 
the unphysical region biases the pdf and produces a sys- 
tematic error. The remedy to this is the method proposed 
by Feldman & Cousins (1998). We, however, cannot apply 
their method to the data we use, it has to be done at the 
time of original data analysis. We shall only give mention 
where it has been done, and where it should have been 
done. 

Since none of the likelihood contours in the data we 
use look normal nor even symmetric, we clearly choose 
the maximum likelihood method. (The tool to use is ac- 
tually not the likelihood function which is the product of 
individual pdf 's, but the negative of the sum of their log- 
arithms, or the log-likelihood function.) We approximate 
the contours by general fifth order polynomials of the form 

P(Q ro Q A ) = QM, m + n<5. (2) 
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There are then 20 terms in the polynomials, so we read 
off 20 points from the la, 1.64(7, 2a, 3a contours and 
the best value, where available. Since we already know 
the approximate location of the globally favored region 
from all the previous studies, it is enough that we require 
our polynomial approximation to be good over that re- 
gion. This fit region is defined by 0.15 < TO < 0.50 and 
0.40 < Oa < 0.88, but the sample points are of course 
taken also from outside this region in order to obtain a 
well-behaved polynomial inside the region. Far away from 
it the polynomial approximation of course breaks down 
completely. 

One should recognize that the observational likelihood 
surfaces in the Q m ,Q,\ plane are not known with a very 
good resolution. Thus one cannot set very high require- 
ments on the polynomial representation inside the la con- 
tours. We have been checking that the polynomial is non- 
negative in the fit region, and that its minimum (i.e. of 
the negative log-likelihood function) is indeed at the ob- 
servational best value, where reported. We have also been 
checking the location of the 0.33ct contour of the polyno- 
mial approximation, in order to verify that it is reasonably 
centrally located with respect to the la contour and to the 
best value, when reported. 

3. DATA 

Altogether we use six independent data sets meeting our 
criteria, grouped into SN la data, CMB data, LSS data 
and other data. But as we shall see, some of these data sets 
actually comprise several other important observations. 

3.1. SUPERNOVA la DATA 

The SN la observations by the High-z Supernova Search 
Team (HSST) of Riess et al. (1998) and the Supernova 
Cosmology Project (SCP) of Perlmutter et al. (1998, 1999) 
are well enough known not to require a detailed presen- 
tation here. The importance of these observations lies in 
that they determine approximately the linear combination 
Q\ — fl m which is orthogonal to Qo = + ^A- 

HSST use two quite distinct methods of light-curve fit- 
ting to determine the distance moduli of their 16 SNe la 
under study. Their luminosity distances are used to place 
constraints on six cosmological parameters: h, Q m , CIa, qo, 
and the dynamical age of the Universe, to- The MLCS 
method involves statistical methods at a more refined level 
than the more empirical template model. The moduli are 
found from a x 2 analysis using an empirical model con- 
taining four free parameters. The MLCS method and the 
template method give moduli which differ by about la. 
Once the distance moduli are known, the parameters h, 
Sl m , fl\ are determined by a maximum likelihood fit, and 
finally the Hubble constant is integrated out. (The results 
are really independent of h.) One may perhaps be some- 
what concerned about the assumption that each distance 



modulus is normally distributed. We have no reason to 
doubt that, but if the iterative \ 2 analysis has yielded sys- 
tematically skewed pdf 's, then the maximum likelihood fit 
will amplify the skewness. 

The authors state that "the dominant source of sta- 
tistical uncertainty is the extinction measurement" . The 
main doubt raised about the SN la observations is the 
risk that (part of) the reddening of the SNe la could be 
caused by intervening dust rather than by the cosmolog- 
ical expansion. Among the possible systematic errors in- 
vestigated is also that associated with extinction. No sys- 
tematic error is found to be important here, but for such 
a small sample of SNe la one can expect that the selection 
bias might be the largest problem. 

The authors do not express any view about which 
method should be considered more reliable, thus noting 
that " we must consider the difference between the cosmo- 
logical constraints reached from the two fitting methods 
to be a systematic uncertainty". We shall come back to 
this question later. Here we would like to point out that 
if one corrects for the unphysical region fl m < using the 
method of Feldman & Cousins (1998), the best value and 
the confidence contours will be shifted slightly towards 
higher values of Q m . This shift will be more important for 
the MLCS method than for the template method, because 
the former extends deeper into the unphysical region. 

The likelihood contours in the fl m , ft a plane (Riess et 
al. 1998, Fig. 6) correspond to 68.3%, 95.4% and 99.7% 
confidence, respectively (la, 2a, 3a). 

Let us now turn to SCP, which studied 42 SNe la. The 
MLCS method described above is basically repeated, but 
modified in many details for which we refer the reader to 
the source. The distance moduli arc again found from a 
X 2 analysis using an empirical model containing four free 
parameters, but this model is slightly different from the 
HSST treatment. The parameters fl m and J7a are then de- 
termined by a maximum likelihood fit to four parameters, 
of which the parameters M.$ ( an absolute magnitude) and 
a (the slope of the width-luminosity relation) are just an- 
cillary variables which are integrated out (h does not en- 
ter at all). The authors then correct the resulting likeli- 
hood contours for the unphysical region TO < using the 
method of Feldman & Cousins (1998). The likelihood con- 
tours in the Cl m , J7a plane (Perlmutter et al. 1999, Fig. 7) 
correspond to 68% (la), 90% (1.64a), 95% (1.96a), and 
99% (2.58cr) confidence, respectively. 

Since the number of SNe la is here so much larger than 
in HSST, the effects of selection and of possible systematic 
errors can be investigated more thoroughly. SCP quotes 
a total possible systematic uncertainty to £^ a * and Q^ at 
of 0.05. 

If we compare the observations along the line defining 
a flat Universe, SCP finds fl A - fl m = 0.44 ± 0.085 ± 0.05, 
whereas HSST finds Q A - fl m = 0.36 ± 0.10 for the 
MLCS method and fi A - M m = 0.68 ± 0.09 for the tem- 
plate method. This comparison tells us that the template 
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method is afflicted by systematic errors of its own. We 
choose the former since this method is basically the same 
as used by SCP. This is admittedly a selection bias of 
ours, but we shall account for it, at least partly, by apply- 
ing the same systematic error of Af^ at = Af2 A ia * — 0.05 
to HSST as to SCP. 

SCP and HSST then agree within their statistical er- 
rors - how well they agree cannot be established since 
they are not completely independent. Part of the differ- 
ence may be explained by selection bias in the smaller set 
of SNe la. 

In Fig. 1 we show the confidence contours of the log- 
likelihood sum of the two SN la observations in our poly- 
nomial approximation, drawn only in the ranges of Q m 
and f^A that we sample. Along the flat line these obser- 
vations determine Qa — = 0.45 ± 0.13. Note that our 
value is not obtained from the weighted mean of the SCP 
and HSST values, but from our log-likelihood sum. 

There are many other types of observations which give 
complementary information in support of the SN la data. 
These observations have been summarized briefly by Perl- 
mutter et al. (1998, 1999), and in more detail for instance 
by us (Roos & Harun-or-Rashid 1998, 1999, 2000). But 
there have also been gravitational lensing data in strong 
conflict with the SN la data. The best value in our Fig. 1 is 
excluded with 99.7% confidence by the joint optical (spi- 
ral galaxies) and radio data of Falco, Kochanek & Muoz 
(1998) (six gravitational lenses analyzed). However, the au- 
thors point out that these results depend on the choice of 
galaxy sub-type luminosity functions in the lens models. 
Subsequently Chiba & Yoshii (1999) have emphasized this 
point in an analysis with E/S0 luminosity functions that 
yielded a best fit mass density in a flat cosmology, finding 
^A - = 0.4 + 0.2/ - 0.4 in agreement with the SN la 
data. 

More recently Helbig (2000) has shown preliminary 
results from the Cosmic Lens All-Sky Survey (CLASS) 
of radio lenses only. These results are still inaccurate, 
—0.8 < fl\ — < 0.3 at 95% confidence, but appear 
to be in strong conflict with the SN la data. It is still too 
early to say whether SNe la will have to come down to 
smaller values of Oa — H m , or whether gravitational lens- 
ing will have to go up. In Section 5 we shall discuss how 
our combined fit changes when including the preliminary 
CLASS constraints. 

3.2. CMB DATA 

Before the advent of the balloon observations BOOMER- 
ANG (de Bernardis et al. 2000) and MAXIMA-1 (Balbi 
et al. 2000, Hanany etal. 2000), Linewcavcr (1998) and 
Tegmark (1999) analyzed all the then existing CMB data 
in the form of multipole power spectra up to I ~ 800. 
The parameter space is then very large, but some pa- 
rameters really drop out and others can be handled in 
a simplified manner if their effect is significant only be- 
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Fig. 1. The confidence contours of the log- likelihood sum of 
the two SN la observations (HSST and SCP). The curves cor- 
respond to lcr and 2<j in the (r2 m , fiA)-plane. The significance 
of the square is described in the text. The diagonal line corre- 
sponds to a flat cosmology. 



low or above t ~ 100. In the analysis of Tegmark (1999), 
the following ten cosmological parameters are jointly con- 
strained: n k ,fl^, the optical depth parameter r, the am- 
plitudes and slopes A s , n s ,A t ,n t of scalar and tensor fluc- 
tuations, and the physical matter densities Ub , ^ c dm , oj v . 
Of these parameters only six are well constrained; the 
resulting 6-dimensional likelihood function is then inte- 
grated over remaining parameters, and it is stated to be 
highly non-Gaussian in some directions. We trust that the 
plotted marginalized 2-dimensional confidence limits in 
the O to ,Oa plane are then realistic, and do not contain 
any imposed Gaussian form. 

The approximations made are claimed to reproduce 
the power spectrum to about 5% accuracy. Otherwise no 
systematic errors are discussed. But the input data show 
rather large scatter, so one might hope that by combining 
them, most of the systematic differences between them 
have been taken into account. One worry is that the best 
fitting models all fail to quite match the COBE DMR data 
(Tegmark 1999). 

The likelihood contours in the fl m , f^A plane from this 
compilation (Tegmark 1999, Fig. 3) correspond to 68% 
(lcr) and 95% (1.96<r) confidence, respectively. No best 
value point is given. 

The balloon observations BOOMERANG (de Bernardis 
et al. 2000) and MAXIMA-1 (Balbi et al. 2000, Hanany 
et al. 2000) have produced the first high-resolution, 
high signal-to-noise maps of the CMB from independent 
patches of the sky, and thereby derived the angular power 
spectrum with a better precision than was achievable in 
compilations of earlier observations. In both cases the so 
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far published results represent a complete analysis of a 
limited portion of the data. 

Given the multipole spectrum, Balbi et al. (2000) fit 
different 7-dimcnsional CDM models to some pixelization 
of the measured range of /, including the COBE DMR 
data. The jointly constrained cosmological parameters are 
t, fih, fi m , Oa, n s and Cio, the amplitude of fluctuations at 
multipole I = 10. For a seventh parameter they alterna- 
tively used h and flbh 2 . Marginalizing over five parame- 
ters, the confidence range in the fi m , ft a plane were found. 
The first acoustic peak is the dominant feature in the 
power spectrum, the maximum occurring at £ = 197 ±6 in 
BOOMERANG and at I ~ 220 in MAXIMA- 1. The po- 
sition of the peak determines Oo pretty independently of 
all other parameters; thus for a flat universe it determines 
f2 TO . 

BOOMERANG (de Bernardis et al. 2000) jointly 
constrain a six-dimensional parameter space, comprising 
h, flbh 2 , fi m , Q\, n s and an overall normalization A. Sub- 
sequently they marginalize over four parameters to obtain 
the confidence range in the fl m ,fl\ plane plotted in their 
Fig. 3. 

The difference in n between the two BOOMERANG 
and MAXIMA-1 is rather large, since I ~ 200ftj7 1-58 for 
TO = 0.3 and fio near 1 (Weinberg 2000). (In the litera- 
ture one sometimes sees the relation I ~ 200Sl f 7°' 5 which 
is true only when CI a = 0.) This difference clearly repre- 
sents a systematic error which should be allowed to affect 
the total fit. All other systematic errors discussed have 
less influence on £l m and ilo- The likelihood contours in 
the fi m ,fi A plane from MAXIMA-1 (Balbi et al. 2000) 
correspond to 68% (la), 95% (1.96cr), and 99% (2.58er) 
confidence, respectively, which we can use. No best value 
point is given. 

BOOMERANG has only published a coarsely pixelized 
likelihood surface of 95% (1.96c) confidence (de Bernardis 
et al. 2000, Fig. 3), which would make our polynomial 
fit a poor approximation. Therefore we do not include 
BOOMERANG as an independent constraint. However, 
it is included together with the LSS constraint to be dis- 
cussed in the next subsection, so the BOOMERANG in- 
formation is not neglected. 

A totally different approach is taken by Jaffe et al. 
(2000). The two data sets both have a calibration uncer- 
tainty, 20% for BOOMERANG and 8% for MAXIMA-1, 
which Jaffe et al. (2000) uses to adjust the power spec- 
tra so that the peaks are more similar in amplitude. The 
data can then be combined into multipole bands, and the 
goodness-of-fit improves considerably. However, since the 
relative importance of the two peaks is then altered, the 
combined data yields a significantly shifted Sl m , and the 
originally visible systematic difference in Q TO disappears. 
We fear that this leads to an underestimation of the Q, m 
systematic error. As stressed before, we prefer to exhibit 
all systematic errors in order to achieve minimum loss of 
information. 



3.3. LSS AND OTHER DATA 

An important source of information on TO and flo is the 
power spectrum of matter density fluctuations. Durrer & 
Novosyadlyj (1999) have chosen to study this for Abell- 
ACO clusters (Retzlaff et al., 1998) , arguing that the 
power spectrum of clusters should better represent the 
observed Universe as a whole, rather than density fluctu- 
ations on the scale of galaxies. Making this choice may 
represent a bias implying ignoring some inherent system- 
atic error, but this seems the only choice for us at the 
moment. 

Novosyadlyj et al. (2000a) combine the power spec- 
trum of Abell-ACO clusters with six independent con- 
straints for the amplitude of the fluctuation power spec- 
trum on different scales: from clusters at different red- 
shifts, from quasar spectra, and from the bulk flow of 
galaxies in our vicinity (cf. their reference list). Very im- 
portantly, they also use the CMB value I = 197 ± 6 
for the multipole moment of the first acoustic peak from 
BOOMERANG which thus gets included into our data 
base. In addition they constrain the Hubble constant to 
be h = 0.65±0.10, a low but reasonable compromise value 
with a generous error. They constrain the baryon density 
to be fl b h 2 = 0.020 ±0.002 (95%CL) (Buries et al. 2001). 

The MAXIMA-1 and BOOMERANG values for Q b h 2 
are in stark contradiction to that value, but this conflict 
does not affect the value of Slo noticeably. In fact, the out- 
come of these authors' multiparameter least squares fit to 
the cluster power spectrum yields a baryon density more 
similar to the MAXIMA-1 and BOOMERANG value, and 
yields a higher Hubble constant closer to present best de- 
terminations. 

It is not clear whether any systematic errors are in- 
cluded by Novosyadlyj et al. (2000a). However, the fit 
errors are rather large anyway, e.g. TO = 0.37+qis 
marginalized over all other parameters, so that the in- 
fluence of systematic errors would be small. The informa- 
tion in Novosyadlyj et al. (2000b), Fig. la are the la, 2a 
and 3a confidence contours in the fi m , f2 A plane, obtained 
by marginalization over the other five parameters (and 
assuming that only one species of massive neutrinos con- 
tribute to the neutrino density parameter fl v ). 

Our last constraint is based on the measurements of 
the coordinate distance to sources at redshifts between 1 
and 2, which depends on the global values of the cosmo- 
logical parameters, but which is independent of the power 
spectrum of density fluctuations, and of whether matter 
is biased relative to light. From a parent population of 70 
powerful extended classical double radio galaxies Guerra 
et al. (2000) have studied a subset of 20, for which it was 
possible to estimate independently the mean and charac- 
teristic size, and thus the coordinate distance. There is 
one model parameter (3 in the theory in addition to Q m 
and Oa- The main systematic error is the model uncer- 
tainty in (3, but it is shown to be unimportant compared 
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to known statistical errors, and quite uncorrelated to TO 
and Oa- 

The information in Guerra et al. (2000), Fig 11, are 
the 68% (la) and 90% (1.64ct) likelihood contours in the 
Q m ,Q\ plane, obtained by marginalizing over j3. The fa- 
vored region is quite large, so that this constraint is at 
present quite weak, adding only some preference for small 
tt m values. A best value is claimed at Q m = and 
Oa = 0.45. It is not clear to us how the unphysical re- 
gion fl m < has been treated, in any case not with the 
Fcldman - Cousins (1998) procedure. 

In Fig. 2 we show the confidence contours of the log- 
likelihood sum of the data sets discussed in this subsection 
and in the CMB subsection. We have plotted our polyno- 
mial approximation only in the ranges of f2 m and that 
we sample. As can be clearly seen, the likelihood func- 
tion contains information mainly on fio, but it also gives 
a rather conspicuous upper limit on the orthogonal com- 
bination Oa — ^m- 




Fig. 2. The confidence contours of the log-likelihood sum 
of MAXIMA-1, early CMB, Double Radio Galaxies, and 
LSS (including BOOMERANG and various other indepen- 
dent constraints). The curves correspond to la and 2a in the 
(n m , f2A)-plane. The significance of the square is described in 
the text. The diagonal line corresponds to a flat cosmology. 



4. RESULTS 

Adding up our polynomial approximations to the confi- 
dence contours of all the data discussed in the previous 
section, results in Fig. 3, where we show the location of 
the minimum, the la and 2a contours. From this Figure 
one can read off the following results: 



0.31 



+0.12 
-0.09 



n A = 0.68 ±0.12, 



(3) 
(4) 



or alternatively 

Q = 0.99 ±0.04 



n A - n m = 0.37 



+0.20 
-0.23- 



(5) 
(6) 



Of these results, only the determination of fig is quite 
precise and worth detailed attention. We can conclude 
from it that a flat universe with Sl = 1 is very likely. 

To the results in Fig. 3 we have to add some further 
quantifiable systematic errors which we evaluate as fol- 
lows. 

As mentioned earlier, Perlmutter et al. (1998, 1999) 
have quoted a total systematic error for n^ at and fl f A lat 
along the flat line of ±0.05. We consider that the same 
error should be applied to the SN la data of HSST, where 
a similar evaluation did not give a significant result due 
to the limited sample of SNe la. As explained in Sec. 3.1, 
this is further motivated by the discord between the MLCS 
method and the template method. Displacing both the SN 
la contours by ±0.05 along the flat line, we obtain a very 
small systematic error in the fio direction 

Ai(O )= ±8;SS1- (7) 

There are also two kinds of systematic errors inher- 
ent to our method of analysis. Firstly we are reading off 
the coordinates of the confidence contours of the different 
observations with some finite precision. We estimate this 
precision to be 

A 2 (n ) = 0.027. (8) 

Secondly, since we only use 20 points to fit the confi- 
dence contours of each observations, there is an arbitrari- 
ness in their choice; all we require is that the confidence 
contours should be well fitted by whichever polynomial. 
We have tested this polynomial arbitrariness and found 
that it results in the systematic error 



A 3 (O ) = 0.01. 



(9) 



The quadratic sum of the errors in Eqs. (7), (8), (9) is 
then 

A tot (O ) - 0.03. (10) 

Thus our final result for Q is 

ft = 0.99 ±0.04 ±0.03, (11) 

where the first error is statistical and the second error 
systematical. Thus our total error is ±0.05. 

Let us now turn to the case of exact flatness, fl m = 
1 — Q\. Along the flat line the SN la systematic error is 

Ai(fi ) /,ot = ±0.025. (12) 

Our result is then 

n^ at = 0.31 ±0.04 ±0.04, (13) 

where the first error is statistical and the second error 
systematical. Thus our total error here is ±0.055. Note 
once again that this systematic error is not included in 
Fig.3. 
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5. DISCUSSION 

The closest comparison we can make with other anal- 
yses of Oo is that of Jaffe et al.(2000), who combine 
BOOMERANG, MAXIMA-1, COBE DMR, the SN la 
constraints of SCP (but not HSST), and another LSS in- 
put than we have used. Also, the earlier CMB data com- 
piled by Lineweaver (1998) and Tegmark (1999), and the 
information from the double radio sources are not used. 
They find 

Q = 1.06 ±0.04. (14) 

to be compared with our result in Eq. (11). The difference 
in central value is easy to understand, and due to two 
causes. The main cause is the way BOOMERANG and 
MAXIMA-1 have been combined, as we explained already 
in the Data section. A small shift in the same direction is 
due to our inclusion of other input. The statistical errors 
are the same, but we have in addition a systematic error, 
part of which is not applicable to the Jaffe et al.(2000) 
analysis. 

A rather similar analysis is that of Durrer & Novosyad- 
lyj (2000), who find n m + fl A « 1.06 and fl k w -0.06, 
where our definition of firj corresponds to Cl m + ft a + &k = 
1. On the flat line these authors find fl m = 0.35 ± 0.05. 
We also note that these authors have taken into account 
all the observational constraints we used in our previous 
analyses (Roos & Harun-or-Rashid 1998, 1999, 2000), and 
which we therefore did not refer to explicitly here. 

One further constraint which we did not make use of 
here is the position of the peak in the matter power spec- 
trum of quasars at z « 2 as observed by Roukema & Ma- 
mon (2000). The reason for the omission is that their like- 
lihood contours in the fi m , Qa space are so jagged that our 
polynomial approximations just cannot reproduce them. 
To quote one result, they find 

n^ at = 0.30 ±0.15 (15) 

in good agreement with us. The error here is so large that 
the inclusion of this result would not have changed our 
conclusions. 

As we mentioned in Sec. 3.1, Helbig (2000) (fig. 3) 
has plotted preliminary lensing constraints from the Cos- 
mic Lens All-Sky Survey (CLASS) of radio lenses, which 
appear to be in strong conflict with the SNe la data. A 
joint analysis is statistically quite meaningless, but ex- 
cluding one or the other is a biased choice. Since we have 
so far included the SNe la data and excluded the lens- 
ing data, let us now include both. The result is shown 
in Fig. 4. One notes that the best value then moves 
to n rn = 0.34+°;^, fl A = 0.63±° %, or alternatively 
O = 0.97t°;^, fl A - fl m = 0.29±°;^. We note that this 
best value is excluded by the SNe la data at ltr CL, and 
it is excluded by the lensing data at 97% CL. 

We conclude from Sec. 4 that tto is equal to unity to 
within ±0.05 and that ftf at = 0.31 to within ±0.055. 




0.2 0.4 0.6 0.8 1 



Fig. 3. Figs. 1 and 2 combined. The '+' marks the best fit: 
(£2 m , ^a) = (0.31,0.68). The diagonal line corresponds to a flat 
cosmology. 
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Observation 


Reference 


Source 


SN la: HSST 


Riess et al. (1998) 


Fig. 6 


SN la: SCP 


Perlmutter et al. (1999) 


Fig. 7 


CMB: MAXIMA-1 


Balbi et al. (2000) 


Fig. 3 


CMB compilation 


Tegmark (1999) 


Fig. 3 


Double Radio Gal. 


Guerra et al. (2000) 


Fig. 11 


LSS 


Novosyadlyj et al. (2000b) 


Fig. 1(a) 


LENSING 


Helbig (2000) 


Fig. 3 



Table 1. The observational data used in the fifth order poly- 
nomials are summarized. 
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